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, In the framework of a microscopic string model inclusive charged parti- 

■ cle distribution and baryon and antibaryon production are described. The 

C*") ' emphasis is put on high energies (RHIC) where shadowing corrections play 

a crucial role. Some recent developments on J ftp suppression at CERN-SPS 
are also discussed. Possible consequences for the crucial issue of thermal 
equilibration of the produced system are considered. 



1. Introduction 



This work is a continuation of the one presented in ref pQ , where I dis- 
cussed multiparticle production at CERN-SPS energies. The present article 
is mainly concerned with higher energies where the efffects of shadowing 
play a very important role. In the framework of the model presented be- 
low, the shadowing corrections can be computed from high mass diffraction 
practically without any new free parameter. When the effects of shadowing 
are taken into account, the model describes the inclusive charged particle 
production at RHIC as a function of centrality 2 . A comparison with 
the results obtained [3J 0] in the framework of the saturation model is also 
presented. 

Another new development concerns net baryon production (stopping). 
It is shown that both SPS and RHIC data can be described with the same 
mechanism (and the same values of the parameters) used in pp interactions. 
This indicates that there is no evidence for an "anomalous" stopping in the 
heavy ion data jSJ. 

As already observed at CERN-SPS, rare processes like strange and mul- 
tistrange baryon and antibaryon production, can only be described with the 
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introduction of some final state interaction between the produced particles 
(comovers interaction) 1 j. It turns out, however, that the interaction cross- 
sections required to describe the data are comparatively small (a few tenths 
of a mb) and, in view of the shortness of the interaction time (5 -j- 7 fm) it 
seem quite improbable that the system can reach thermal equilibrium. In a 
recent development [5] reported below, we show that the same formalism of 
final state interaction used at CERN-SPS can describe RHIC data with the 
same values of the parameters. Predictions for S and E production have 
been confirmed by recent STAR data. Predictions for £1 and f2 are also 
given. 

Finally, we analyze the new NA50 data on J/tfj suppression at CERN- 
SPS in the comovers approach [7] [S] and discuss expectations at RHIC. 

2. The model 

2.1. Hadron-Hadron Interactions 

The Dual Parton Model (DPM) [9 and the Quark Gluon String Model 
(QGSM) TO] are closely related dynamical models of soft hadronic interac- 
tions. They are based on the large-iV expansion of non-perturbative QCD 
[11-13] and on Gribov's Reggeon Field Theory (RFT) |14j . Their main aim 
is to determine the mechanism of multiparticle production in hadronic and 
nuclear interactions. The basic mechanism is well known in e + e~ annihila- 
tion (Fig. 1). Here the e + e~ converts into a virtual photon, which decays 
into a qq pair. In the rest system of the virtual photon the quark (with 
colour 3) and the antiquark (colour 3) separate from each other producing 
one string (or chain) of hadrons, i.e. two back-to-back jets. Processes of 
this type are called one-string processes. 

In hadron-hadron interactions, a one-string mechanism is also possible 
but only in some cases, namely when the projectile contains an antiquark 
(quark) of the same type than a quark (antiquark) of the target, which 
can annihilate with each other in their interaction. For instance in n + p, 
the d of 7r + can annihilate with the d of p and a single string is stretched 
between the u of ir + (colour 3) and a diquark uu of p (colour 3). This 
mechanism is also possible in pp interactions (Fig. 2) but not in pp. This 
already indicates that it cannot give the dominant contribution at high 
energy. Indeed, when taking the square of the diagram of Fig. 2 (in the sense 
of unitarity) we obtain a planar graph, which is the dominant one according 
to the large- ./V expansion. However, this only means that this graph has the 
strongest coupling. Since flavour quantum numbers are exchanged between 
projectile and target, this graph gives a contribution to the total cross- 
section that decreases as an inverse power of s {1/yfs). A decrease with s is 
always associated with flavor exchange. For instance, the charge exchange 
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7r~p — ► 7r°n cross-section also decreases as \j\fs. Actually, the diagram in 
Fig. 2 corresponds to the exchange of a secondary Reggeon, with intercept 
close to 1/2. 

In order to prevent the exchange of flavour between projectile and tar- 
get, the d and d have to stay, respectively, in the projectile and target 
hemispheres. Since they are coloured, they must hadronize stretching a sec- 
ond string of type d-d. We obtain in this way a two-string diagram (Figs. 
3-5)). Taking the square of this diagram, we obtain a graph with the topol- 
ogy of a cylinder (Fig. 6). It turns out that this is the simplest topology 
one can construct which does not vanish as s —* oo due to flavour exchange. 
Therefore, we obtain in this way the dominant graph for hadron-hadron 
scattering at high energy. The diagram in Fig. 6 is called a Pomeron (P) 
and the graphs in Figs. 3-5 a cut Pomeron. Its order in the large- N expan- 
sion is 1/N 2 [12-13]. Note that due to energy conservation the longitudinal 
momentum fractions taken by the two systems at the string ends have to 
add up to unity. 

There are also higher order diagrams (in the sense of the large- ./V ex- 
pansion) with 4, 6, 8 strings which give non-vanishing contributions at high 
energy. An example of the next-to-leading graph for pp interactions is shown 
in Fig. 7. It contains four strings - the two extra strings are stretched be- 
tween sea quarks and antiquarks. The square of this graph corresponds to 
a PP cut and has the topology of a cylinder with a handle. Its order in the 
large- N expansion is l/N 4 . The one with six strings corresponds to a PPP 
cut and to the topology of a cylinder with two handles (order 1/N e ), etc. 

The single particle inclusive spectrum is then given by [S] 



dN PP 



(y) = E — E °n (N^ (y) + N*>-**(y) + (2n - 2)iV« (y 



dy 

~ N™-«°(y) + JV*-«(y) + (2k - 2)N^(y) (1) 
where k = J2 na n / a n is the average number of inelastic collisions. Note 

n n 

that each term consists of 2n strings, i.e. two strings per inelastic collisions. 
Two of these strings, of type qq-q, contain the diquarks of the colliding 
protons. All other strings are of type q-q. 

The weights a n of the different graphs, i.e. their contribution to the total 
cross-section, cannot be computed in the large-A^ expansion. However, it 
has been shown |15j that there is a one-to-one correspondence between the 
various graphs in the large- N expansion and those in perturbative Reggeon 
Field Theory |14j . We use the weights obtained from the latter - with the 
parameters determined from a fit to total and elastic cross-sections jH] , ^0] . 
At SPS energies we get k = 1.4 and at RHIC k = 2 at y/s = 130 GeV and 
k = 2.2 at yfe = 200 GeV 0. 
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The hadronic spectra of the individual strings N(y) are obtained from 
convolutions of momentum distribution functions, giving the probability 
to find a given constituent (valence quark, sea quark of diquark) in the 
projectile or in the target, with the corresponding fragmentation functions. 
The dependence of N(y) on the number of collisions appears via the former. 
It is a result of energy conservation. (The larger the number of strings, the 
smaller the average invariant mass of each one). 

Momentum distribution and fragmentation functions are largely deter- 
mined from known Regge intercepts jU] jlUj . The momentum distribution 
function of a valence quark in a hadron behaves as 1/y/x. As in the parton 
model, this behaviour results from the intercept 1/2 of a Reggeon trajectory. 
Thus, in average, the valence quark in a proton is slow and the diquark is 
fast due to energy conservation. Both momentum distribution and fragmen- 
tation functions are assumed to be universal, i.e. the same in all hadronic 
and nuclear interactions. This property gives to the model a great predic- 
tive power. Finally, individual strings are assumed to be independent. In 
this way, the hadronic spectra of a given graph are obtained by adding up 
the corresponding ones for the individual strings. This leads to a picture, 
in which, for any individual graph, particles are produced with only short- 
range (in rapidity) correlations. Long-range correlations (and a broadening 
of the multiplicity distributions) are due to fluctuations in the number of 
strings, i.e. to the superposition of different graphs with their correspond- 
ing weights. This gives a simple and successful description of the data in 
hadron-hadron and hadron- nucleus interactions 0| ^Oj |16| . 



2.2. Nucleus-Nucleus Interactions 

The generalization of Eq. to nucleus-nucleus collisions is rather 
straighforward. For simplicity let us consider the case of A A collisions and 
let n A and n be the average number of participants of each nucleus and 
the average number of binary NN collisions, respectively. At fixed impact 
parameter b, we have ^7] 

^(6) = n A (b) [NfifM + N^\y) + (2k - 2)N^ 
+ {n(b)-n A {b))2kN q ;-^{y) (2) 

n A (b) and n(b) are computed from the standard formulae in the Glauber 
model. The physical meaning of Eq. @ is quite obvious. The expression in 
brackets corresponds to a NN collision. Since n A nucleons of each nucleus 
participate in the collision, this expression has to be multiplied by n A . Note 
that in Eq. Q the average number of collisions is k — and the number 
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of strings 2k. In the present case the total average number of collisions 
is kn - and the number of strings 2 /ere. The second term in Eq. is 
precisely needed in order to have the total number of strings required by 
the model. Note that there are In a strings involving the valence quarks and 
diquarks of the participating nucleons. The remaining strings are necessarily 
stretched between sea quarks and antiquarks. The value of fj,(b) is given 
by fx(b) = kv(b) with v(b) = n(b)/riA(b) (i(b) represents the total average 
number of inelastic collisions suffered by each nucleon (for more details see 
Section 5). 

We see from eq. (j2J) that dN AA /dy is obtained as a linear combination of 
the average number of participants and of binary collisions. The coefficients 
are determined within the model and depend on the impact parameter via 
n(b). Note that the presence of a term proportional to the number of bi- 
nary collisions is a general feature of RFT and is not related to minijet 
production. 

As discussed in Section 2.1 the average invariant mass of a string con- 
taining a diquark at one end is larger than the one of a q-q string since 
the average momentum fraction taken by a diquark is larger than that of 
quark. It turns out that the same is true for the central plateau, i.e. : 
N qq ~ q {y* ~ 0) > N q ~ q (y* ~ 0). Let us now consider two limiting cases : 

dN AA 

If N qs - q °(y* ~ 0) < N qq - qv (y* ~ 0) , then — — (y* ~ 0) ~ n A ~ A 1 (3) 



If N^-^iy* ~ 0) ~ N qq - q Hy* ~ 0) , then (y* ~ 0) - n ~ A 4 / 3 . 

ay 

(4) 

In the first case we obtain a proportionality in the number of participants 
tla whereas in the second case we obtain a proportionality in the number of 
binary collisions. Since dN AA /dy = (1 / 'a AA)do AA / 'dy , the latter result im- 
plies that da AA /dy ~ A 2 , i.e. all unitarity corrections cancel and we obtain 
the same result as in the impulse approximation (Born term only). This re- 
sult is known as the Abramovsky-Gribov-Kancheli (AGK) cancellation and 
is valid for a general class of models which includes the Glauber and eikonal 
ones. It implies that, for the inclusive cross-section, soft and hard processes 
have the same A-dependence. However, the AGK cancellation is violated 
by diagrams related to the diffraction production of large-mass states - the 
so-called triple Pomeron or enhanced diagrams. These diagrams give rise 
to shadowing corrections as discussed below. Their effect is very important 
in nuclear collisions since they are enhanced by A 1 / 3 factors. 



G 



2.3. Shadowing corrections 

In Appendix A, we discuss the physical content of the AGK cutting 
rules and their practical realization in the probabilistic Glauber model. It is 
shown there that multiple scattering diagrams, resulting from the s-channel 
iteration of the Born term, give non-vanishing contributions to the total 
cross-sections (shadowing). However, in the case of the single particle inclu- 
sive cross-section, these contributions cancel identically (AGK cancellation), 
provided the measured particle has been produced in an inelastic interac- 
tion (cut Pomeron). If, on the contrary, the trigger particle is produced 
in the vertex function (blob) of the multiple scattering diagram, one ob- 
tains the same shadowing effects than in the total cross-section. This is the 
physical origin of the AGK violations present in DPM (see section 2.2). It 
is clear that if the blob has a small extension in rapidity, production from 
the blob will mainly contribute to the fragmentation region. Therefore, 
at mid-rapidities, and sufficient large energy, the AGK cancellation will be 
valid. 

Let us consider next the contribution to the total cross-section resulting 
from the diffractive production of large mass states. Clearly, this is equiv- 
alent to an increase of the rapidity extension of the blob - which, in this 
case, can cover the mid-rapidity region. Therefore, shadowing corrections 
to the single particle cross-section will be present in this case, provided the 
measured particle is part of the diffractively produced system. As shown in 
Appendix A, the shadowing correction is just given by the diffractive cross- 
section with negative sign. (This is exactly true only for purely imaginary 
amplitudes). The theoretical expression of the diffractive cross-section is 
well-known. An important part is given by the triple-Pomeron term. It has 
also been measured experimentally and, thus, the shadowing corrections can 
be computed with no free parameters. 

Considering for simplicity only the contribution of the triple Pomeron 
term, the effect of the shadowing corrections is obtained |181 12] by multi- 
plying Eq. (j2J) by 

p ,,x fd 2 s f A (s) f B (b-s) 

RAB(b) = « (5) 

where 

Here the function F(s) is given by the integral of the ratio of the triple 
Pomeron cross-section d?a ppp /dydt at t = to the single Pomeron ex- 
change cross-section a p (s) : 
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rVmax 

F{s) = 4vr / dy 

" Xlm.im 



1 ffo 



(Tp(s) dy dt 



C [exp (Ay ) - exp (Ay min )} 



t=o 



(7) 

with y = £n(s/M 2 ), where M 2 is the squared mass of the diffractive sys- 
tem. For a particle produced at y cm = 0, y m ax = \&n{s / mr) 2 and y m in = 
en{R A m N /V3). A = a P (0) - 1 = 0.13 and C is a constant proportional 
to the triple Pomeron coupling. Ra is the nuclear radius, 7a (6) the nuclear 
profile function and Tab(P) = J sTa(s)Tb(6 — s). 

Eqs. to l|7 |) can be derived only when the triple Pomeron coupling is 
small and, thus, the second term in the denominator of © is small compared 
to the first one. In this case, we have [1 + AF(s)Ta(o)] ~ 1 — AF(s)Ta(P), 
and only the contribution of the triple Pomeron graph is involved in the 
shadowing. In general higher order rescatterings are also present. They are 
model dependent. The denominator in Eq. © correspond to the sum of all 
"fan" diagrams with Pomeron branchings (generalized Schwimmer model 

B9D- 

It is interesting to study the A-dependence of the shadowing corrections 
in the limit of large triple Pomeron coupling (when the first term in the 
denominator of © can be neglected). In this case we find Raa ~ A~ 2 ^ s , 
i.e. the A 4 / 3 behaviour resulting from the AGK cancellation is reduced to 
A 2 / 3 . This limit was considered by Kancheli many years ago [2U] . 

Note that shadowing corrections to inclusive spectra are not especific to 
soft processes. The triple Pomeron terms described above are also respon- 
sible for shadowing in hard processes. 



3. Charged Particle Multiplicities 

3.1. Low pt 

At SPS energies the limit given by Eq. @ is not reached, and Eq. © 
leads to an A dependence of dN AA /dy at y* ~ in A a with a only slightly 
above unity, (a ~ 1.08 between 2 and 370 participants). On the other hand, 
shadowing corrections are small due to phase Space limitations (jjmax ~ Drain 
in Eq. Q). The results g] for Pb Pb collisions at ^/i = 17.3 GeV are 
shown in Fig. 8. We see that both the absolute values and the centrality 
dependence are well reproduced. When the energy increases, Eq. © shows 
that the value of a should increase towards 4/3, in the absence of shadowing 
corrections. However, the effect of the latter is increasingly important and, 
as a result, the value of a varies little with s. At y/s = 130 GeV, without 
shadowing corrections the A-dependence is A a , with a ~ 1.27 in the same 
range of n par t - a value which is not far from the maximal one, a = 4/3 from 
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Eq. (j3J). With the shadowing corrections the A-dependence is much weaker 
(lower line of the shaded area in Fig. 9) [2J. The ^4-dependence is now given 
by A a with a ~ 1.13 - always in the range of n par t from 2 to 370. As we see, 
the increase of a from SPS to RHIC energies is rather small. This value of 
a is predicted to change very little between RHIC and LHC, where a ~ 1.1. 
For, the increase from a ~ 1.27 to a ~ 4/3 obtained in the absence of 
shadowing is compensated by an increase in the strength of the shadowing 
corrections, leaving the effective value of a practically unchanged. 

3. 2. Large pr 

Let us define the ratio 

_ <PN AA /dy dp T 
KAA(b) ~ n{b) d?NPP/dy dp T ■ (8j 

For central Au Au collisions at y* ~ 0, we find R AA = ^ 113 ~ 4 / 3 = 0.34 
when numerator and denominator are integrated over pr. Clearly this value 
corresponds to small values of pr which give the dominant contribution to 
dN/dy. This result is in agreement with the measured values of Raa at 
Pt ~< Pt > pT] . This was to be expected from the results presented in 
Fig. 9. It is interesting that these data, as well as PHENIX ones |22j show 
approximately the same value of Raa at large pr- 1 

More precisely, the data show a small increase of Raa at moderate pr 
and, at large pr, they show a scaling in the number of participants - rather 
than in the number of binary collisions. Such a result is expected in the 
present approach. Indeed, at large pr the shadowing corrections strongly 
decrease due to the increase of mr in y max (eq. ©)• However, the larger 
threshold at large pr affects mostly the q-q strings - which have a smaller 
invariant mass than the qq-q ones. Thus, with increasing pr we approach 
the limit in eq. (jSJ) leading to a scaling in the number of participants. The 
small increase of Raa at moderate pr is probably due to the Cronin effect 
or to a combination of this effect and jet quenching. Recently, it has been 
shown [22] that the Cronin effect is considerably smaller at RHIC and LHC 
energies than at SPS ones due to the change with energy of the coherence 
length. 

4. Comparison with the saturation model 

In the saturation model, the ^-dependence of charged particle inclusive 
spectrum in the saturation regime (i.e. Aqcd "C Pt < Qs, where Q s is the 



1 I would like to thank N. Armesto, K. Boreskov, Y. Dokshitzer, A. Kaidalov, O. 
Kancheli, A. Krzywicki and D. Schiff for discussions on this subject. 
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saturation scale) is given by [31 E] 

dN ^ _^/__ ^ 

Comparison with previous section results indicates that, apart from loga- 
rithmic factors due to a s (Q s ), we obtain the same result as in RFT with 
MAXIMAL shadowing. As discussed in the previous section, this result is 
in violent disagreement with RHIC data. 

The question is then how a reasonable description of the data has been 
obtained in The answer is the following. The authors have considered 
dN/dy rather than dN / dyd 2 px- By integrating over d 2 pT up to Q s and 
assuming a p^-broadening corresponding to Q 2 ~ A 1 /^ they have gained 
one power of A 1 ^. Furthermore, the factor aj 1 ~ log A 1 /^ increases the 
effective power slightly above A 1 reaching agreement with experiment. The 
problem with this explanation is that a p^-broadening in A 1 / 3 is much 
larger than the one seen in the data (which is of the order of 30 % between 
peripheral and central collisions). In order to describe it a parametrization 
of the saturation scale has been introduced in [3] of the type : 



0.61 + 0.39 



^ n par t(b) V 3 



Tlpart- 



max 



Pso ■ (10) 



With this definition, p so is the value of the saturation scale for the most 
central collisions (corresponding to n par t-max = 347 participants). With 
such an expression of the saturation scale, the ^-dependence of ^ is the 

same as that of dN/dyd 2 px within 30 % and, with the mild ^4-dependence 
of a" 1 used before, it is not possible to describe the data. In view of that 
the authors use instead the following expression 



a s 1 ~ log 



0.61 + 0.39 




(11) 



A 2 

where fi 2 = Q ° D is chosen to be 0.6. Note the fine tuning between the 0.61 

Ps,0 

in the numerator and 0.6 in the denominator of (|11|) . As a consequence, the 
value of a^ 1 f° r peripheral collisions is extremely small and a^ 1 increases 
by a large factor between peripheral and central collisions. With this fine 
tuning agreement with experiment is recovered. Note that for peripheral 
collisions (n par t(b) ^ n vart - m ax) one is practically sitting on the Landau 
pole (i.e. the argument of the log is very close to unity). Note also that 
with Aqcd = 200 MeV, the value of the saturation scale for the most central 
collisions is very small (260 MeV). 
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In the saturation model the ^4-dependence of dN/dyd 2 pT grows larger 
with increasing pj- and it turns out that the model can reproduce [21j the 
measured values of the ratio Raa, Eq. ©, in the range 1.5 GeV < pr < 
5 GeV. However, as discussed above the model has an ^4-dependence at 
lower pt which is too weak. Thus, it is perhaps not surprising that it has 
the right ^4-dependence for larger values of pt in some range. 

The above considerations indicate that saturation is not reached at 
RHIC energies. The considerations in the previous section based on RFT 
suggest that it will not be reached at LHC either. 

5. Nuclear Stopping Revisited 

In pp collisions the net proton (p-p) distribution is large in the fragmen- 
tation regions and has a deep minimum at mid-rapidities. In contrast to this 
situation a much flatter distribution has been observed |22] in central Pb 
Pb collisions at CERN-SPS 2 . In view of that, several authors have claimed 
that the stopping in heavy ion collisions is anomalous, in the sense that it 
cannot be reproduced with the same mechanism (and the same values of 
the parameters) used to describe the pp data. In a recent paper [S] it has 
been shown that this claim is not correct. 

In the model described in previous sections, the net baryon can be pro- 
duced directly from the fragmentation of the diquark. Another possibility is 
that the diquark splits producing a leading meson in the first string break- 
up and the net baryon is produced in a further break-up. Clearly, in the 
first case, the net baryon distribution will be more concentrated in the 
fragmentation region than in the second case. The corresponding rapidity 
distributions are related to the intercepts of the relevant Regge trajectories, 
a qq and a q , respectively, i.e. they are given by e Ay ^ a \ Here Ay is the 
difference between the rapidity of the produced net baryon and the maxi- 
mal one. In the case of the first component, in order to slow down the net 
baryon it is necessary to slow down a diquark. The corresponding Regge 
trajectory is called baryonium and its intercept is known experimentally to 
be a qq = —1.5 ± 0.5. For the second component, where a valence quark is 
slowed down, we take a q = 1/2 3 . 

In this way we arrive to the following two component model for net 
baryon production out of a single nucleon 

2 Actually, a huge stopping was first observed at AGS. However, in this case we are in 
a different regime (intra- nuclear cascade). 

3 There is a third possibility in which the net-baryon transfer in rapidity takes place 
without valence quarks (string junction or gluonic mechanism) with intercept either 
as j = 1/2 26 or otsj = 1 |27|. We find no evidence for such a component from the 
existing pp and AA data. Its smallness could be related to the fact that it produces 
an extra string of hadrons and, thus, does not correspond to the dominant topology. 
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M (y) = aC Z 1-a9(0) (l - ^ + )/^W-3/2+n„K(o)-^(o)) 
dy M + 

+ (1 - a)C; Z+~^ (0) x (1 - Z+ )M6)~3/2+ c +n s ,K(0)-«,(0) (12) 

where n sq is the number of strange quarks in the hyperon a p (0) = 1/2 
0(f)(0) = 0, Z + = (e y ~ ymax ) ? y max is the maximal value of the baryon rapidity 
and fi(b) is the average number of inelastic collisions suffered by the nucleon 
at fixed impact parameter b (see section 2.2). The constants C p and 
are obtained from the normalization to unity of each term. The small Z 
behaviour is contro lied by the corresponding intercept. The factor (1 — 
£ + ) M (6)-3/2 ig obtained by requiring that the Z-fractions of all quarks at 
the ends of the strings, other than the one in which the baryon is produced, 
go to zero [9-10]. Following conventional Regge rules an extra a p (0) = 
a^(0) = 1/2 is added to the power of 1 — Z + for each strange quark in the 
hyperon. 

The fraction, a, of the DB breaking component is treated as a free 
parameter. The same for the parameter c in the DP component - which 
has to be determined from the shape of the (non-diffractive) proton inclusive 
cross-section in the baryon fragmentation region. It can be seen from Eq. 
(|12|) that stopping increases with /i(6), i.e. with the total number of inelastic 
collisions suffered by each nucleon. This effect is present in the two terms of 
(|12l) and is a consequence of energy conservation. The question is whether 
or not this "normal" stopping is sufficient to reproduce the data. In other 
words whether or not the data can be described with a universal value of a, 
i.e. independent of fi and the same for all reactions. 

Eq. ()12|) gives the total net baryon density, but it does not allow to 
determine the relative densities of different baryon species. In order to do 
so we use the simple quark counting rules described in Appendix B. 

A good description of the data on the rapidity distribution of pp — > 
p-p + X both at -v/s = 17.2 GeV and y/s = 27 A GeV is obtained from Eq. 
(|12|) with a = 0.4, c = 1, a q = 1/2 and a gq = —1. The results are shown 
in Table 1 at three different energies, and compared with the data. As we 
see the agreement is reasonable. For comparison with the nucleus-nucleus 
results, all values in Table 1 have been scaled by the number of participants 
pairs in central Pb Pb and Au Au collisions (n^ = 175). As it is well known, 
a pronounced minimum is present at y* = 0. There is also a substantial 
decrease of the mid-rapidity yields with increasing energy. Also, the mid- 
rapidity distributions get flatter with increasing energy since the net proton 
peaks are shifted towards the fragmentation regions. 

It is now possible to compute the corresponding net baryon production 
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in heavy ion collisions and to check whether or not the data can be described 
with Eq. (|12[) using the same set of parameters as in pp. 

The results [3] for net proton (p — p) and net baryon {B — B) production 
in central Pb Pb collisions at y^s = 17.2 GeV and central Au Au collisions at 
y/s = 130 GeV are given in Table 2. The centrality is defined by the average 
number of participants - n par t = 2ua = 350 in both cases. Experimental 
results are given in brackets. 

The comparison of column 2 with the pp results in Table 1 at the same 
energy, shows the well known change in the shape of the rapidity distribution 
between pp and central Pb Pb collisions at SPS. The minimum at y* = 
is much less pronounced in Pb Pb and the net proton peaks in the pp 
fragmentation regions are shifted to y* ~ ±1.5. More interesting are the 
results in columns 4 and 5 which contain the predictions for Au Au at RHIC. 
We see that the shape of the rapidity distribution is very different from the 
one at SPS. 

In conclusion, we have found that "anomalous" stopping is not needed 
in order to describe the present data. Related models have been proposed 
in [29-31]. The results for heavy ion collisions are rather similar to the ones 
obtained from eq. 1)12(1 . However, in these models there is some increase in 
the size of the second component with the number of inelastic collisions, i.e. 
some anomalous stopping is present. 

6. Hyperon-Antihyperon Production 

Strange particle production and, in particular, of multistrange hyperons 
has been proposed as a signal of Quark Gluon Plasma formation. Flavor 
equilibration is very efficient in a plasma due to large gluon densities and 
low thresholds. An analysis of the results at SPS in the framework of the 
present model has been presented in [3T]. In the following we concentrate 
on RHIC results. 

A general result in DPM is that the ratios B /h~ and B /h~ of baryon 
and antibaryon yields over negatives decrease with increasing centralities. 
This is easy to see from Eq. (J2J. The production from q s — q s strings 
scales with the number of binary collisions. These strings have a smaller 
(average) invariant mass than the qq-q strings and, thus, are more affected 
by the thresholds needed for BB pair production. As a consequence, the 
centrality dependence of B and B production will be smaller than the one 
of negatives. (The same effect was discussed in Section 3.2 in connection 
with large pt production). The effect is rather small at RHIC energies. 
However, it is sizable and increases with the mass of the produced baryon. 
In contrast with this situation, the data for A's show no such decrease and 
an increase is present for 3 production. Data on production are not yet 
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available. However, SPS data clearly show a hierachy in the sense that the 
enhancement of baryon production increase with the mass (or strange quark 
content) of the produced baryon. 

The only way out we have found is to give up the assumption of string 
independence. Until now we have assumed that particles produced in dif- 
ferent strings are independent from each other. In the following we allow 
for some final state interactions between comoving hadrons or partons (see 
Section 8). We proceed as follows. 

The hadronic densities obtained in Section 2 are used as initial con- 
ditions in the gain and loss differential equations which govern final state 
interactions. In the conventional derivation [22] of these equations, one uses 
cylindrical space-time variables and assumes boost invariance. Furthermore, 
one assumes that the dilution in time of the densities is only due to longitu- 
dinal motion 4 , which leads to a r _1 dependence on the longitudinal proper 
time r. These equations can be written [312] [3J 



The first term in the r.h.s. of (|13j) describes the production (gain) of particles 
of type i resulting from the interaction of particles k and t. The second term 
describes the loss of particles of type i due to its interactions with particles 
of type k. In Eq. (|TH|) pi = dNi/dyd 2 s(y,b) are the particles yields per 
unit rapidity and per unit of transverse area, at fixed impact parameter. 
They can be obtained from the rapidity densities © using the geometry, 
i.e. the s-dependence of ha and n. The procedure is explained in detail in 
[7] where the pion fragmentation functions are also given. Those of kaons 
and baryons can be found in UJ. ctm are the corresponding cross-sections 
averaged over the momentum distribution of the colliding particles. 

Equations ()13|) have to be integrated from initial time tq to freeze-out 
time tj. They are invariant under the change r — > ct and, thus, the re- 
sult depends only on the ratio tj/tq. We use the inverse proportionality 
between proper time and densities and put tj/to = (dN/dyd 2 s(b))/pf. 
Here the numerator is given by the DPM particles densities. We take 
pj = [3/TrRp](dN~ /dy) y *^Q = 2 fm~ 2 , which corresponds to the density of 
charged and neutrals per unit rapidity in a pp collisions at yfs = 130 GeV. 
This density is about 70 % larger than at SPS energies. Since the corre- 
sponding increase in the AA density is comparable, the average duration 

4 Transverse expansion is neglected. The fact that HBT radii are similar at SPS and 
RHIC and of the order of magnitude of the nuclear radii, seems to indicate that this 
expansion is not large. The effect of a small transverse expansion can presumably be 
taken into account by a small change of the final state interactions cross-sections. 




(13) 
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time of the interaction will be approximately the same at CERN SPS and 
RHIC - about 5 to 7 fm. 

Next, we specify the channels that have been taken into account in our 
calculations. They are 

ttN ^ KA(T,) , vrA(S) ^ KE , irE ^ KVL (14) 
We have also taken into account the strangeness exchange reactions 

vrA(S) KN , ttE ^ KA(S) , ntl ^ KE (15) 

as well as the channels corresponding to (|14|) and Q15|) for antiparticles 5 . 
We have taken = a = 0.2 mb, i.e. a single value for all reactions in (|14[) 
and (JT5|) - the same value used in ref. [2j to describe the CERN SPS data. 

Before discussing the numerical results and the comparison with exper- 
iment let us examine the qualitative effects of comovers interaction. As 
explained in the beginning of this Section, without final state interactions 
all ratios K/h~, B/h~ and B/h~ decrease with increasing centrality. The 
final state interactions (|14|). ()15|) lead to a gain of strange particle yields. 
The reason for this is the following. In the first direct reaction (|14j) we have 
pw > PKi Pn > PA, PnPN 3> PkPA- The same is true for all direct reaction 
(|14|) . In view of that, the effect of the inverse reactions (|14j) is small. On 
the contrary, in all reactions (|15|) . the product of densities in the initial 
and final state are comparable and the direct and inverse reactions tend to 
compensate with each other. Baryons with the largest strange quark con- 
tent, which find themselves at the end of the chain of direct reactions (j!4j) 
and have the smallest yield before final state interaction, have the largest 
enhancement. Moreover, the gain in the yield of strange baryons is larger 
than the one of antibaryons since ps > pg- Furthermore, the enhancement 
of all baryon species increases with centrality, since the gain, resulting from 
the first term in Eq. (|13|) . contains a product of densities and thus, increases 
quadratically with increasing centrality. 

6.1. Numerical Results 

All our results refer to mid-rapidities. The calculations have been per- 
formed in the interval —0.35 < y* < 0.35. In Fig. lOa-lOd we show the ra- 

5 To be precise, of all possible charge combinations in reactions 1141 . we have only kept 
those involving the annihilation of a light q-q pair and production of an s-s in the 
s-channel. The other reactions, involving three quarks in the f-channel intermediate 
state, have substantially smaller cross-sections and have been neglected. All channels 
involving 7r° have been taken with cross-section a/2 since only one of the uu and dd 
components of 7r° can participate to a given charge combination. For details see the 
first paper of |31|. 
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pidity densities of B, B and B — B e versus h~ = dN~ /drj = (I /1. 17) dN/dy 
and compare them with available data [33-35]. We would like to stress 
that the results for H and 3 were given [B] before the data This is an 
important success of our approach. 

In first approximation, the yields of p, p, A and A yields over h~ are 
independent of centrality. Quantitatively, there is a slight decrease with 
centrality oipjhr and pjh~ ratios, a slight increase of A/h~ and A/h~ and 
a much larger increase for 3 (H)//i _ and O (£l)/h~. This is better seen 
in Fig. 11a and lib where we plot the yields of B and B per participant 
normalized to the same ratio for peripheral collisions versus n par t- The 
enhancement of B and B increases with the number of strange quarks in 
the baryon. This increase is comparable to the one found at SPS between pA 
and central Pb Pb collisions, especially for antibaryons. The ratio K~/tc~ 
increases by 30 % in the same centrality range, between 0.11 and 0.14 in 
agreement with present data. The ratios B/B have a mild decrease with 
centrality of about 15 % for all baryon species - which is also seen in the 
data. Our values for N ch /N^ ax = 1/2 are : p/p = 0.69, A/A = 0.74, 
3/3 = 0.79, O/Q = .83, to be compared with the measured values [3B] : 



p/p = 0.63±0.02±0.06 , A/A = 0.73±0.03 , 3/3 = 0.83±0.03±0.05 . 

The ratio K + /K~ = 1.1 and has a mild increase with centrality, a feature 
also seen in the data. 

Note that a single parameter has been adjusted in order to determine the 
absolute yields of BB pair production, namely the p one - which has been 
adjusted to the experimental p value for peripheral collisions. The yields of 
all other BB pairs has been determined using the quark counting rules given 
in Appendix B. The experimental data in Fig. 10 are not corrected for feed- 
down from weak decays. If these corrections were the same (in percentage) 
for all baryon species, our results should be compared with uncorrected 
yields. This seems to be the case for p, p, A and A where the feed-down 
corrections are of the order of 20 %. As a consequence, our predictions for 
3, 3, and have a 20 % uncertainty. 

Although the inverse slopes ("temperature") have not been discussed 
here, let us note that in DPM they are approximately the same for all 
baryons and antibaryons both before and after final state interaction - the 
effect of final state interaction on these slopes being rather small |37| . 



In the numerical calculations the net baryon yields have been obtained using the 
approach in |S| and |31l . This approach is conceptually different from the one in 
Section 5 but the numerical results are similar. 
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7. New J/ip suppression data and the comovers interpretation 

The NA38-NA50 collaboration have observed a decrease of the ratio of 
J/ip to dimuon (DY) cross-sections with increasing centrality in SU and 
Pb Pb collisions. The same phenomenon has been observed in pA collisions 
with increasing values of A. In this case, it is interpreted as due to the 
interaction of the pre-resonant cc pair with the nucleons of the nucleus it 
meets in its path (nuclear absorption). As a result of this interaction, the 
cc pair is modified in such a way that, after interaction, it has no projection 
into J/ip (a DD pair is produced instead). The J/tp survival probability 
S a bs is well known (see for instance eq. (7) of |7j) and depends on a single 
free parameter a a b s , i.e. the absorptive cc — N cross section. 

The NA50 collaboration has shown that the J/ip suppression in Pb Pb 
collisions has an anomalous component, i.e. it cannot be reproduced using 
nuclear absorption alone. Two main interpretations have been proposed : 
deconfinement and comovers interaction. The latter mechanism has been 
described in Section 6 for strange particle production. In the case of J/ip 
suppression, a single channel is important namely cc (or J/ip) interacting 
with comoving hadrons and producing a DD pair. In this case, eq. (|13|) 
can be solved analytically. The expression of the survival probability S co 
can be found in (see eq. (8)). It depends on a free parameter a co , i.e. 
the effective cross-section for the comovers interaction. 

Two important sets of new data have been presented recently by the 
NA50 collaboration on pA and Pb Pb collisions jSH] • Before these data 
were available, the NA50 interpretation of the data was as follows. The pA, 
SU and peripheral Pb Pb data can be described with nuclear absorption 
alone, with a a bs = 6.4 ± 0.8 mb. At Et ~ 40 GeV there is a sudden onset 
of anomalous suppression with a steady fall off at large Et- However, at 
variance with this view, the most peripheral Pb Pb points lied above the 
nuclear absorption curve - which extrapolates pA and SU data. 

The new pA data indicate a substantially smaller value of the absorptive 
cross-section. However, within errors, pA and SU data can still be described 
with a a bs = 4.4 ± 0.5 mb |38| . The new Pb Pb preliminary data, taken in 
2000 with a target under vacuum, are consistent with previous ones except 
for the most peripheral ones - which are now lower and consistent with the 
nuclear absorption curve JSH] • in this way, the NA50 interpretation remains 
valid. However, the new data lend support to the interpretation based on 
comovers interaction. Indeed, due to the smaller value of a a b s there is more 
room for comovers interaction (i.e. for anomalous suppression) in SU. 

Actually, before the new data were available, it has been argued |S] that 
a value of a a bs = 4.5 mb is also consistent within errors with the old pA data. 
Using this value and a co = 1 mb it has been possible to describe all available 
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data within the comovers scenario 003 ED]- There was, however, a caveat, 
as pointed out in Indeed, there was a mismatch of about 30 % between 
the absolute normalizations in SU and Pb Pb. Actually, the ratio of the first 
normalization to the second one is only 1.04 ± 0.02 |38| . (This factor takes 
into account both the isospin correction in SU and the rescaling in energy). 
This mismatch was induced by the high values of the most peripheral Pb 
Pb data in the old NA50 data. Indeed, since the relative contribution of the 
comovers to J/ip suppression is larger for central collisions, the centrality 
dependence of the J/ip suppression gets flatter (steeper) with decreasing 
(increasing) values of a co (at fixed cr a fe s ). In order to reproduce the shape 
of the J /if) over DY ratio, from very peripheral to central collisions, in 
the old NA50 analysis, a value of <r co = 1 mb was required. On the other 
hand the new data are described with a smaller value a co = 0.65 mb. This 
decrease of a co leads to a decrease of the absolute normalization, which is 
now consistent with the SU one 7 . 

The results |42) of the comover interaction model with a a b s = 4.5 mb 
and a co = 0.65 mb are presented in Fig. 12. As in ref. [3U] the steady 
fall-off of the J/ip over DY ratio at large Et is obtained introducing the 
Et fluctuations. The agreement with the new NA50 data is quite 
satisfactory. The absolute normalization is 47. The corresponding one in 
SU is 45 in perfect agreement with the expectations discussed above. 

It is interesting to note that the data obtained using the Et calorimeter 
and the zero-degree calorimeter (ZDC) analysis are consistent with each 
other when using the measured Et — Ezdc correlation. This result was 
predicted in ref. jSJ. 

Next, I would like to discuss briefly the expectations for J/ip supression 
at RHIC in the comovers interaction model. The calculation of the survival 
probability S co is quite safe. Indeed, since a co is a cross-section near thresh- 
old, the same value obtained at SPS should be used at RHIC. The situation 
is quite different for S a bs- Many authors assume that a a f, s is the same at 
RHIC and at SPS. It has also been suggested that it can be significantly 
larger at RHIC. However, it seems plausible that at mid-rapidities, nuclear 
absorption at RHIC is small due to the fact that, contrary to SPS, the cc 
pair is produced outside the colliding nuclei. It is therefore crucial to have 
data on J /if) production in pA interactions at RHIC. If S a b s ~ 1 the J/ip 
suppression at RHIC and SPS will be comparable since the smallness of the 
nuclear absorption will be approximately compensated by the increase of 
the comovers suppression - due to a larger comovers density at RHIC. Very 



It is interesting that almost the same value of o co (p co = 0.7 mb) was obtained in 
|41| from an analysis of SU data and old Pb Pb data (which covered a much smaller 
centrality range). In |41| the absolute normalization in SU and Pb Pb were in good 
agreement with each other. 
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preliminary data tend to indicate that this is indeed the case. 

A quantitative analysis of the new NA50 data in the deconfining sce- 
nario is still missing. On the other hand, the centrality dependence of the 
average px of J/tp is better described in the comovers approach than in a 
deconfining scenario [13]. At RHIC energies, a small nuclear absorption in 
pA collisions (i.e. S a b s ~ 1), would be a very interesting situation in order 
to discriminate the comovers interaction model from a deconfining scenario. 
Indeed, in the latter, the shape of the centrality dependence would be al- 
most flat for peripheral collisions (below the deconfining threshold) and 
would decrease above the threshold. Such a behaviour would be a clear sig- 
nal of deconfinement. On the contrary, in the comovers scenario, the fall-off 
would be continuous, from peripheral to central collisions, and determined 
by the same value of a co obtained from CERN SPS data. 

8. Conclusions 

Quark Gluon Plasma (QGP) formation is obtained in statistical QCD, 
i.e. QCD applied to a system in thermal equilibrium. Therefore, one of 
the main issues in heavy ions physics is to determine whether or not the 
produced final state reaches thermal equilibrium. An argument in favor of 
equilibrium is the fact that particle abundances are well described in terms 
of statistical models. However, one should take into account that statistical 
models are also very successful in pp and even e + e~ interactions. Therefore, 
it is important to study whether or not particle abundances can be obtained 
in a microscopic model such as DPM. 

As a starting point we have assumed that particles produced in different 
strings are independent (see Section 2). In this case thermal equilibrium 
cannot be reached no matter how large the energy density is. Indeed, in 
this case a large energy-density is the result of piling up a large number 
of independent strings. The assumption of independence of strings works 
remarkably well in hh and hA interactions jH] ^U] - even in the case of 
event samples with 5 or 6 times the average multiplicity - indicating that 
no sizable final state interaction is present in these reactions. In nucleus- 
nucleus collisions, we have described charged particle inclusive production 
and its centrality dependence. The model exhibits a term proportional to 
the number of binary collisions - which has been seen in the data both at 
SPS and RHIC. The presence of such a term is required by unitarity and is 
not due to minijets. 

However, it is clear that in heavy ion collisions, where several strings 
occupy a transverse area of 1 fm 2 , the assumption of string independence 
has to break down. This is indeed the case. As we have seen, some data 
cannot be described without final state interaction. It could have happened 
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that this final state interaction is so strong that the string picture breaks 
down and becomes totally useless. This does not seem to be the case. 
On the contrary, present data can be described using the particle densities 
computed in the model as initial conditions in the gain and loss (transport) 
equations governing the final state interaction. The interaction cross-section 
turns out to be small (of the order of a few tenths of a mb). Due to this 
smallness and to the limited interaction time available, final state interaction 
has an important effect only on rare processes, in particular H, O and J/ip 
production. The bulk of the final state is not affected. 

Of course it is not possible to conclude that thermal equilibrium has not 
reached. However, particle abundances not only do not allow to conclude 
that it has been reached, but, on the contrary, their centrality dependence 
tends to indicate that this is not the case. Let us consider for instance p 
and p production. In our model, their yields are practically not affected by 
final state interaction, i.e. they are practically the same assuming string 
independence. Yet, the model reproduces the data, from very peripheral 
to very central interaction. This success would be difficult to understand 
in a QGP scenario in which for peripheral collisions (below the critical 
density) there is strong, non-equilibrated, pp annihilation, which becomes 
equilibrated for central ones, above the critical density. More generally, the 
QGP scenario would be strongly supported if some kind of threshold would 
be found in the strange baryon yields around the critical density value. 
At SPS energies, evidence for such a threshold in the S yield has been 
claimed by the NA57 collaboration 44 . Unfortunately, these data only 
cover a limited range of centrality. In contrast to this situation the RHIC 
data explore the whole centrality range from very peripheral to very central 
collisions and the centrality dependence of the yields of p, A, 5 and their 
antiparticles shows no structure whatsoever. If the same happens for £1 and 
production (as predicted in our approach) the case for QGP formation 
from strange baryon enhancement will be rather weak. 

Finally, it should be stressed that the final state interaction of comovers 
in our approach is by no means a trivial hadronic effect. Indeed, the inter- 
action of comovers starts at the early times where densities, as computed in 
DPM, are very large. In this situation the comovers are not hadrons (there 
are several of them in the volume occupied by one hadron, and, moreover, 
at these early times hadrons are not yet formed). This is probably the rea- 
son why in our approach the comover interaction cross-sections required to 
describe the data are smaller than in a hadron gas model. 

Before concluding, I would like to say that it is an honor and a pleasure 
to contribute to this special issue of Acta Physica Polonica in hommage to 
my friend Jan Kwiecinski. I would also like to acknowledge his important 
contributions to the model presented here, realized during his (too rare) 
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visits to Orsay. In particular he played an important role in the general- 
ization of the Dual Parton Model to heavy ions collisions |17j and also in 
introducing 03] a semi-hard component in the model. 

Appendix A 



a) Reggeon Field Theory versus Glauber Model 

The reggeon calculus or reggeon field theory (RFT) ^1] provides a field 
theoretical formulation of the eikonal (for hh collisions) or the Glauber (for 
hA and AB) models, valid at high energies. The main difference between the 
RFT and the Glauber model is that, at high energies, the coherence length 
is large and the whole nucleus is involved in the interaction. Moreover, due 
to the space-time development of the interaction, when, at high energy a 
projectile interacts inelastically with a nucleon of the nucleus, the formation 
time of (most of) the produced particles is larger than the nuclear size and, 
thus, particles are produced outside the nucleus. Therefore, planar diagrams 
give a vanishing contribution at high energy. The relevant diagrams are non- 
planar, describing the "parallel" interactions of constituents of the projectile 
with the target nucleons (in the case of an hA collision). This picture is in 
clear contrast with the Glauber model, in which the projectile undergoes 
successive (billiard ball type of) collisions with the nucleons of the target. 

In spite of these differences, one recovers the Glauber formula in first 
approximation. This formula corresponds to the contribution of the initial 
state (on-shell projectile pole) to the various rescattering terms. In RFT 
one has, besides these contributions, also the contributions due to low mass 
and high mass diffractive excitations of the projectile. The latter are very 
important since, as we have seen in Section 2.3, they give rise to shadowing 
corrections. 

b) Cutting Rules 

An important feature of RFT is that it obeys to the so-called AGK 
cutting rules |46| . These rules allow to relate to each other the different s- 
channel discontinuities of a given graph, and also to relate them to the con- 
tribution of this graph to the total cross-section. In this way, they provide 
a powerful link between total cross-section and multiparticle production. In 
order to illustrate these rules, let us consider the case of an interaction of 
a hadron h with two different nucleons of the target nucleus A (with A — 2 
spectators), and let us assume that the object exchanged in the i-channel 
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of each collision is purely imaginary (Pomeron). 

Let us consider the cutting by a plane in between the two interactions 
(i.e. in between the two Pomerons). We obtain in this way a diffractive 
intermediate state containing a large rapidity gap. Let us call +1 its contri- 
bution to a tot ■ From the cutting rules we find that the inelastic contribution 
obtained cutting through one of the interactions (an interference term) has 
a weight —2 relative to the previous one. Since there are two interactions 
one can cut through, one obtains —4. Finally, cutting by a plane through 
the two interactions (which is possible since the graph is non-planar) has 
a relative weight +2. This last contribution is also inelastic and has an 
average multiplicity which is twice that of the previous one. The total con- 
tribution of this double scattering, to a tot is thus equal to +1 — 4 + 2 = — 1, 
a negative contribution. The total contribution to the non-diffractive in- 
elastic cross-section is <jnd = — 4 + 2 = —2. We see in this way that the 
(negative) contribution of a double interaction to &nd is two times larger, 
in absolute value that its contribution to a to t ■ In the case of n collisions the 
corresponding factor is 2 n . 

Let us now consider the contribution of a double interaction to the non- 
diffractive single particle inclusive cross-section da/dy. This contribution is 
—4 + 2x2 = 0. Indeed, in the case of the cut through the two interactions 
the contribution to da/dy has an extra factor 2 since the triggered particle 
can be produced in either of the two interactions. It turns out that such a 
cancellation is true to all orders in the number of interactions. We obtain 
in this way the so-called AGK cancellation. All rescattering corrections of 
the Glauber type cancel identically in da/dy. Only the term with a single 
interaction is left - which is proportional to A 1 in pA interactions. 

Note that the crucial ingredient in obtaining the AGK cancellation is 
the fact that the triggered particle has been produced in a cut interaction - 
which gives the extra factor 2 for two cut interactions. The other possibility 
is that, the trigger particle is emitted from the (cut) vertex function (blob). 
Clearly, in this case the extra factor 2 is absent and the AGK cancellation 
is not valid. In this case the shadowing corrections are the same as in the 
total cross-section. 

The AGK cutting rule described above are quite general. They are 
valid in any field theory in which the vertex functions obey the general 
properties of unitarity, crossing and large pr damping. The Glauber model 
is a particular example in which the AGK rules are valid. Their derivation 
in this case is straightforward, as discussed below. 

c) The Probabilistic Glauber Model and the Cutting Rules 

Let us consider for simplicity pA scattering. The main formula of the 
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probabilistic Glauber model is the one that gives the cross-section a n for n 
inelastic collisions of the projectile with n nucleons of the target nucleus, at 
fixed impact parameter b : 

a n (b) = (f) (a inel T A (b)) n (1 - a inel T A {b)) A ~ n (A.l) 

where Oi ne x is the proton- nucleon inelastic cross-section and T A (b) is the 
nuclear profile function. This equation is just the Bernoulli's formula for 
composite probabilities. The first factor is a trivial combinatorial factor 
corresponding to the different ways of choosing n nucleons out of A. The 
second one gives the probability of having n inelastic pN collisions at given 
b. The third one is the probability that the remaining A — n nucleons do 
not interact inelastically. Let us consider first a term with two collisions 
both of which are inelastic. The corresponding cross-section is a 2(b) = 
{2) (cineiTAib)) 2 i.e. a positive term. Let us now consider the case of two 
collisions only one of which is inelastic. The corresponding (interference) 
term is 02(b) obtained from eq. (A.l) by putting n = 1 and taking the 
second term in the expansion of the last factor. We get (r\{b) = —A(A — 
])(&i ne iT A (b)) 2 . We see that a\{b) = —2a\{b). Thus, a rescattering term 
containing two collisions gives a negative contribution to a to t- 

Let us now consider the contribution to da/dy. It is given by a\{b) + 
2a%(b) = 0. Indeed, in the case of a double inelastic collision, the triggered 
particle can be emitted in either of them - hence an extra factor 2. This is 
just the AGK cancellation. It is easy to see that it is valid order by order in 
the total number of collisions. This can also be seen as follows. The total 
inelastic cross-section for pA collision in the Glauber model is given by the 
well known expression 

<nel(b) = E *„(&) = 1 " (1 " Vinel T A {b)) A . (A.2) 

71=1 

This expression contains a term in A 1 (Born term or impulse approxima- 
tion). It also contains contribution from multiple scattering with alternate 
signs. Numerically, it behaves as A a with a ~ 2/3. The single particle 
inclusive cross-section is given by 

daP A A 

— — (b) oc E n Mb) = A a inel T A {b) . (A.3) 

" n=l 

We see that here multiple-scattering contributions cancel identically and 
only the Born term is left. As a consequence of this AGK cancellation 
the ^4-dependence of da/dy in pA interactions behaves as A 1 . In the case 
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of AB collisions it behaves as AB and dN /dy = (\j a AB)da j dy is 
proportional to the number of binary collisions - rather than to the number 
of participants. 

We see in this way that the AGK rules are trivially satisfied in the 
Glauber model. As mentioned in Section A.l in the Glauber model only 
the initial state is present in the vertex function (blob). Thus a secondary 
can only be produced in an interaction and the AGK cancellation is exact. 
In a general theory with a more complicated vertex function, the triggered 
particle may be produced in the blob. As discussed in Section A. 2 this gives 
rise to a violation of the AGK cancellation - which is responsible for the 
shadowing corrections to the inclusive spectra. 
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In order to get the relative densities of each baryon and antibaryon 
species we use simple quark counting rules [SB- Denoting the strangeness 
suppression factor by S/L (with 2L + S = 1), baryons produced out of three 
sea quarks (which is the case for pair production) are given the relative 
weights 

I 3 = 4L 3 : 4L 3 : Y1L 2 S : 3LS 2 : 3LS 2 : S 3 (B.l) 

for p, n, A + S, H°, H - and Q, respectively. The various coefficients of 1% 
are obtained from the power expansion of (2L + S) 3 . 

For net baryon production, we have seen in Section 5 that the baryon can 
contain either one or two sea quarks. The first case corresponds to direct 
diquark fragmentation described by the second term of Eq. (|12|l. The 
second case corresponds to diquark splitting, described by the first term of 
(|12|) . In these two cases, the relative densities of each baryon species are 
respectively given by 

h=L:L:S (B.2) 

for p, n and A + E, and 

I 2 = 2L 2 : 2L 2 : 4LS : -S 2 : -S 2 (B.3) 

2 2 y ' 

for p, n, A + S, H° and S~. The various coefficients in (B.2) and (B.3) are 
obtained from the power expansion of (2L + S) and (2L + S) 2 , respectively. 

In order to take into account the decay of £*(1385) into Air, we redefine 
the relative rate of A's and E's using the empirical rule A = 0.6(E + + E~) 
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- keeping, of course, the total yield of A's plus S's unchanged. In this 
way the normalization constants of all baryon species in pair production 
are determined from one of them. This constant, together with the relative 
normalization of K and tt, are determined from the data for very peripheral 
collisions. In the calculations we use S = 0.1 (S/L = 0.22). 
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Table 1 

Calculated values of the rapidity distribution of pp — > p — p + X at 
v/s = 17.2 GeV and 27.4 GeV (k = 1.4) and y/s = 130 GeV (k = 2). (In 
order to convert da/dy into dN/dy a value of a = 30 mb has been used). 
For comparison with the nucleus- nucleus results, all values in this table have 
been scaled by ha = 175 - the number of participant pairs in central Pb Pb 
and Au Au collisions. Data are in brackets. 



y 


pp — > p — p 
y/s = 17.2 GeV 


PP ~~ * P ~ P 
yfs = 27 A GeV 


pp — > p — p 
= 130 GeV 





9.2 


6.5 


3.6 






[6.3 ± 0.9] 




l 


15.0 


9.3 


4.2 




[16.1 ± 1.8] 


[9.6 ± 0.9] 




1.5 


25.8 


14.6 


5.1 




[24.1 ± 1.4] 


[15.4 ±0.9] 




2 


47.1 


26.2 


6.8 




[45.4 ± 1.4] 


[27.7 ±0.9] 





Table 2 

Calculated values |Sj of the rapidity distribution dN/dy for central Pb 
Pb -^p-p + X and Pb Pb B-B + X &ty^ = 17.2 GeV (k = 1.4) and 
central Au Au -> p-p+V and An An -> B-B+X at ^/s = 130 GeV (fc = 2) 
and y^s = 200 GeV (A; = 2.2) . The centrality has been defined by the 
number of participant pairs (tia = 175 at all energies) and v = n/njy = 4.5, 
5.0 and 5.2 at y/s = 17.2, 130 and 200 GeV, respectively. Data are in 
brackets. 



y 


Pb Pb — > p — p 
y/s = 17.2 GeV 


Pb Pb^ B-B 
y/s = 17.2 GeV 


Au Au — > p — p 
yfs = 130 (200) GeV 


Au Au^B-B 
y/s = 130 (200 GeV) 





23.0 


58.5 


8.0 (7.4) 


20.9 (18.9) 




[26.7 ± 3.7] 


[67.7 ± 7.3] 


[5.6 ± 0.9 ± 24%] 




l 


32.3 


79.7 


9.7 (8.7) 


22.6 (22.0) 




[34.9 ± 1.5] 


[84.7 ± 3.5] 






1.5 


36.3 


87.0 


12.3 (10.9) 


31.5 (27.4) 




[34.4 ± 1.7] 


[80.0 ± 3.9] 






2 


25.3 


57.15 


17.3 (14.3) 


43.4 (35.9) 




[24.7 ± 1.5] 


[56.1 ±3.1] 
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Figure captions 

Figure 1 : The mechanism of particle production in e + e~ annihilation. 
The net of soft gluons and quark loops is only shown here and in Fig. 6. 

Figure 2 : One string diagram in pp. 

Figure 3 : Dominant two-chain (single cut Pomeron) contributions to high 
energy 7r + -proton collisions. 

Figure 4 : Dominant two-chain contribution to proton-antiproton collisions 
at high energies (single cut Pomeron). 

Figure 5 : Dominant two-chain diagram describing multiparticle produc- 
tion in high energy proton-proton collisions (single cut Pomeron). 

Figure 6 : Single Pomeron exchange and its underlying cylindrical topol- 
ogy. This is the dominant contribution to proton-proton elastic scattering 
at high energies. 

Figure 7 : Two cut Pomeron (four-chain) diagram for proton-proton colli- 
sions. 

h igure 8 : The values of dN ch /dy per participant for Pb Pb collisions at 
y/s = 17.3 GeV computed [2j from eq. (2), compared with WA98 data. 

Figure 9 : The values of dN ch /dr] c _ m J (0.5 n par t) for Au Au collisions at 
y/s = 130 GeV computed ;2. from eq. (2) including shadowing corrections 
are given by the dark band in between solid lines. The PHENIX data are 
also shown (black circles and shaded area). 

Figure 10 : (a) Calculated values [5] of dN/dy of p (solid line) p (dashed 
line), and p — p (dotted line) at mid rapidities, |y*| < 0.35, are plotted as a 
function of dN h - /drj, and compared with PHENIX data [33] ; (b) same for 
A and A compared with preliminary STAR data [34] ; (c) same for H~ and 
S compared to preliminary STAR data [35] ; (d) same for 0, and f2. 

Figure 11 : Calculated values 6 of the ratios B/n par t (a) and B/n par t (b), 
normalized to the same ratio for peripheral collisions (n part = 18), plotted 
as a function of n par t. 
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Figure 12 : The ratio of J /ip over DY cross-sections in Pb Pb collisions a 
158 GeV/c versus Et obtained [12] i n the comovers interaction model with 
a abs = 4.5 mb and a co = 0.65 mb. The absolute normalization is 47. The 
preliminary data are from [39]. 
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Figure 6 



Figure 7 




Figure 8 
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Figure 9 
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Figure 10 
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Figure 11a 
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Figure lib 
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